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Abstract. Odd-even effects, also known as staggering effects, are a common feature observed in the yield 
distributions of fragments produced in different types of nuclear reactions. We review old methods, and we 
propose new ones, for a quantitative estimation of these effects as a function of proton or neutron number 
of the reaction products. All methods are compared on the basis of Monte Carlo simulations. We find that 
some are not well suited for the task, the most reliable ones being those based either on a non-linear fit 
with a properly oscillating function or on a third (or fourth) finite difference approach. In any case, high 
statistic is of paramount importance to avoid that spurious structures appear just because of statistical 
fluctuations in the data and of strong correlations among the yields of neighboring fragments. 

PACS. 29.85.Fj Data analysis 


1 Introduction 

An odd-even effect in the yield distributions of final re¬ 
action products has been observed since long time in a 
variety of nuclear reactions. This effect (often called stag¬ 
gering) usually consists in an enhanced production of even 
nuclear species with respect to the neighboring odd ones: 
for example, in the case of charge distributions, even-Z 
species are produced more abundantly than the neighbor¬ 
ing odd-Z one^. 

The staggering effect was first investigated long time 
ago in the fission of actinide nuclei, mainly induced by 
low-energy neutrons (see e.g. [TJ[31[3J|31[3ini[71[5] and refer¬ 
ences therein) and it was attributed to the extra energy 
required to break a pair of protons or neutrons . Experi¬ 
mentally the odd-even effect in fission was found to be less 
pronounced for neutrons than for protons. Many investiga¬ 
tions were devoted to finding a systematic parametrization 
of the normal shapes, against which experimental results 
could be compared to extract the odd-even effects. At very 
low excitation energies, the yield distributions of fission 
fragments are bell shaped, therefore the normal distribu¬ 
tions were assumed to be Gaussians [10]. The odd-even 
effect appeared as a superimposed saw-tooth modulation. 

Having as a reference these semi-empirical shapes, the 
staggering of low-energy fission was usually quantified by 
a single value (called global staggering) extracted from the 
whole fragment distribution. Nowadays staggering effects 
are studied in a variety of different reactions, therefore no 
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^ The opposite effect, i.e. enhancement of odd-Z species and 
reduction of even ones, is rare and is called antistaggering. 


a priori shape of the distribution can be assumed and only 
local values of the staggering can be extracted from small 
adjacent regions of the yield distributions. An extensive 
review of different quantitative methods used for estimat¬ 
ing local and global average values of the odd-even effects 
in fission was given by Gonnenwein m- 

For a quantitative analysis, Amiel and Feldstein m 
US] parametrized the odd-even effects in fragment distri¬ 
butions by assuming that the yields of even-Z fragments 
are intensified by a factor (I -|- A) with respect to the 
smooth distribution YsmiZ) and those of odd-Z fragments 
are reduced by a factor (1 — A), with |A| < 1. The nor¬ 
mal behavior is therefore represented by the arithmetic 
mean between odd and even yield values. The generaliza¬ 
tion, to take into account odd-even effects as a function 
of both proton and neutron numbers (usually assumed to 
be independent from each other), leads to the expression 

YiZ,N) = Y,^{Z,N){l±Az){l±AN) ( 1 ) 

where y^m(Z, N) is the assumed smooth behavior without 
staggering; the sign ‘-I-’ is used for even Z (or N), the 
sign ‘ —’ for odd ones. A = 0 means no odd-even effects, 
while A > 0 indicates enhanced yields of even species and 
reduced yields for odd species (vice versa for A< 0). 

Another approach by Wahl [14] led to a slightly differ¬ 
ent parametrization 

Y{Z,N) = Y,^{Z,N) (Fz)±i (F^)±i , (2) 

with {F)^^ = which can be cast also in the form 

ln(r(Z, N)) = ln(F;^(Z, N)) ±Dz ±Dn (3) 
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where again the sign ‘-I-’ is used for even Z (or N) and 
the sign ‘ —’ for odd ones. In this case F=1 {D=0) corre¬ 
sponds to no odd-even effect, while F>1 {D >0) indicates 
enhanced yields for even species and reduced yields for odd 
species; vice versa for F <1 {D <0). In the parametriza- 
tion by Wahl, the smooth behavior is represented by the 
geometric mean between odd and even yield values. 

Odd-even effects were discovered also in light complex 
fragments produced in high-energy fragmentation and spal¬ 
lation reactions [niiiniiiTiiiHiiiniiiniiii]- The staggering con¬ 
sisted again in an intensification of the yields of even 
charges Z with respect to odd ones, although the reac¬ 
tion mechanism, the shape of the charge distributions and 
the investigated mass region were obviously different from 
those of low-energy fission experiments. 

Recently, staggering effects have been observed also in 
heavy ion collisions from few MeV/nucleon up to Fermi 
energies (15 <E/A< 50 MeV/nucleon) [^1^1^II 25 ll^ 
I271I281I291I30] . These experiments have stirred renewed in¬ 
terest in staggering phenomena. Indeed, in order to study 
the symmetry energy one needs to estimate 

the primary isotopic distributions, which can be reliably 
reconstructed only if the effects of secondary decays are 
small or sufficiently well understood. 

The staggering has been usually ascribed to nuclear 
structure effects that manifest in the decay of the excited 
reaction products or already in the reaction mechanism, 
if part of the reaction proceeds through low excitation 
energies [26]. However, in collisions at intermediate en¬ 
ergies the preferred interpretation is that the odd-even 
staggering effect depends on the structure of the nuclei 
produced near the end of the evaporation chain [191125] . 
The relationship between nuclear structure and odd-even 
yield staggering in nuclear reactions is an active field of 
research. Indeed nice results are emerging in well-chosen 
cases |3S]. However the general case is complex and still 
not completely clarified, as it involves not only the nuclear 
structure of the final nuclei, but also the level densities and 
population probabilities of the parent nuclei. 

Without loss of generality, we will hereafter refer to 
the charge distributions and to the staggering in charge 
Z, but of course exactly the same arguments are valid for 
neutrons too. Often the presence of odd-even staggering 
in Z can be visually appreciated by simple inspection of 
the charge distributions, especially in regions where the 
effect is prominent. For regions with reduced staggering, 
or in the presence of steep variations of the yield as a 
function of charge, the visualization is more difficult and 
some specific handling becomes necessary to highlight the 
existence of staggering effects. 

An objective treatment of the data becomes manda¬ 
tory if one wants to make a quantitative comparison of 
different experiments or different nuclear systems. More¬ 
over, in contrast with fission fragment studies, for other 
reaction mechanisms the shape of the yield distribution 
Ysm has to be deduced from the data with some proce¬ 
dure. The first one, introduced by Tracy et al. |2] in 1972, 
has been widely used for long time (see, e.g., [S]), but 
other treatments have been proposed or can be devised. 


In this paper we examine some methods that could be 
used to put into evidence the presence of the staggering 
and to quantitatively estimate the magnitude of the phe¬ 
nomenon, without relying on an a priori knowledge of the 
yield distribution. These methods are briefly described in 
sect. H] and the simulation in sect. [3] Section 0] presents 
the results and a comparison of the different methods. 

We note that some of the methods presented in this 
paper are commonly used for studying odd-even stagger¬ 
ing on other nuclear quantities, the most notable of all 
being nuclear masses and binding energies (see e.g. |36| 
[37] and references therein). The distinctive feature of frag¬ 
ment yields is that odd-even effects are here a perturba¬ 
tion superposed on a quantity (the yield) that, depending 
on the reaction mechanism, may undergo large and rapid 
changes in magnitude over rather restricted regions of Z. 


2 Smoothing methods 

When analyzing odd-even effects, the main assumption 
of eqs. m-m is that the experimental yield Y{Z) can 
be factorized into the product of the smooth yield Ys^iZ) 
(without staggering effects) multiplied by a staggering fac¬ 
tor. In heavy-ion reaction, where Ysm{Z) is unknown, only 
a local staggering can be determined. The main, essential 
assumption is that the staggering varies with Z in a suffi¬ 
ciently gradual way to be considered - to a good approxi¬ 
mation - constant over an interval of a few Z values. 

The smoothing procedures used for estimating Az or 
Dz from the data are listed below (each with the corre¬ 
sponding abbreviation that will be used in the paper) and 
some of them are schematically illustrated in figs. |T] and 0] 
For more clarity, the procedures that are applied directly 
to the experimental yields will be distinguished from those 
that are applied to their logarithms by prefixing “Y-” in 
the first case and “LY-” in the second case. In both fig¬ 
ures, the open squares joined by dotted lines indicate the 
input cross section without staggering Ysni(Z), while the 
open crosses joined by long dashed lines represent what 
from now on we call the experimental cross section Y(Z), 
because it simulates the outcome of an experiment (for il¬ 
lustration purposes no statistical fluctuations are applied 
in these two figures). The staggering effects are modelled 
d la Amiel and Feldstein with a constant value of 

0.1 for Az- The dotted and long dashed lines are just a 
guide to the eye, joining the points. 

1. Linear interpolation (LIN or Y-2DI). If the shape of 
the charge distribution is sufficiently smooth, the sim¬ 
plest estimate of Ysm{Z) is the average between the 
original experimental distribution and the same dis¬ 
tribution shifted, back and forth, by just one Z unit. 
This procedure tends to eliminate any periodical os¬ 
cillation with period 2Z. Practically, the distribution 
y{Z), which is an estimator of the unknown distribu¬ 
tion Ysui{Z), can be obtained from the experimental 
one with a linear interpolation: 

y{Z) = l{Y{Z-l) + 2YiZ)+Y{Z+l)). (4) 
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In fig. [Ha), y is estimated by means of the solid line 
joining the yields Yz-i/i and Yz+ 1/2 (open diamonds) 
obtained by pairwise averaging the experimental data. 
For each only the measured yields of three consecu¬ 
tive elements (from Z-\ to Z+1) are needed. The stag¬ 
gering parameter jAzj is given by \Y{Z)—y{Z)\ fy{Z) 
and its effect is indicated by the arrows pointing from 
the experimental Y(Z) (crosses) to y(Z) (circles). This 
method is similar to that applied by Zeitlin [TB] (and 
later by others [38l|39]) in high-energy fragmentation 
studies. Note that in his original version, Zeitlin [TB] 
took as a reference for each given Z the line (fine dot¬ 
ted in figure) joining Y(Z—1) and Y(Z+1): in this way 
he obtained an odd-even effect (distance from cross to 
triangle, shown only for Z=7) exactly twice as large as 
the present one (distance from cross to full circle). 

2. Fit-Interpolation over five points (Y-2DIF). In this pa¬ 

per we propose a possible improvement of the pre¬ 
vious method by taking into account also the yields 
of the two neighboring points Y(Z — 2) and Y(Z + 
2). These five yields are pairwise averaged to produce 
four intermediate values Tz- 3 / 2 j T 2 + 1/2 and 

Y 2 + 3/2 (open diamonds), which are used to perform 
a parabolic fit (short dashed line). The value of the fit 
in the central point Z gives the estimated y(Z) and 
the difference with respect to the measured value is 
used to estimate Az- 

3. Parabolic fit (Y-PAR). A new smoothing method was 
recently suggested in Ref. [IH, based on a fitting pro¬ 
cedure too. For each measured value of the yield Y(Z), 
the smoothed value 3 ^(Z) is estimated by fitting a para¬ 
bolic function, 

y(Z) =aZ^ + bZ + c, (5) 

directly to the measured yields over five consecutive 
values of Z, from Z — 2 to Z + 2. Figure [Ijb) illus¬ 
trates this method showing three curves, two for the 
interpolation around even Z (upward parabolas, red 
online) and the other one around an odd Z (down¬ 
ward parabola, blue online). The effect oi Az is again 
indicated by the arrows pointing from the experimen¬ 
tal Y{Z) (crosses) to y(Z) (full circles). The three 
parabolas show the peculiarity of this method. Due 
to the presence of the staggering, the experimental 
points have an up-and-down behavior that is obviously 
not well fitted by a parabola. Therefore the estimated 
effect of staggering (arrows in figure) is enhanced. In 
fact, with respect to the true yields without staggering 
Ysra{Z) (open squares), the estimated values of y(Z) 
(full circles) lie systematically below for even Z values, 
and above for odd ones. This amplifies the staggering 
effect, but does not make the method more sensitive, 
because also the statistical errors are enhanced by the 
same factor. At the same time one has to be aware 
that this amplification hinders a proper quantitative 
comparison with the results of other experiments, eval¬ 
uated with other methods. 



Fig. 1. (color online) Methods Y-2DI and Y-2DIF (a), Y-PAR 
(b) and Y-COS (c). The smooth input distribution (squares) 
is multiplied (crosses) by the staggering factor (1 + ( — 1)^Az)- 
No statistical fluctuations are applied. The full lines represent 
the constructions used to estimate the smoothed yields y{Z) 
(full circles), the arrows show the magnitude of the estimated 
staggering (see text). Panel (a) shows also the method of Zeitlin 
| 18 | (triangle and fine dotted line). 

4. Parabolic fit with oscillations (Y-COS). In this paper 
we propose an improvement of the previous method by 
taking into account the fact that actually, due to stag¬ 
gering, the data oscillate around the assumed smooth 
behavior. The improvement is accomplished by using 
a fit function that consists of a parabola multiplied by 
any properly oscillating function with period 2Z and 
amplitude do, like for example 

Yfit(Z) = (ao Z"^ + bo Z + Co) (1 -l- do cos('kZ)) 

y{Z) (l-bdo cos( 7 rZ)). ( 6 ) 
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Of course any other oscillating function (like a rectan¬ 
gular or a triangular one) with period 2Z and proper 
phase, would have done an equally good job; we choose 
cos(7rZ) because it is simple and easy to implement in 
an analysis code. At the cost of a non-linear fitting pro- 
cedureEl and by adding a fourth parameter, one obtains 
a better and more sensible fit (in fact, Yfit(Z) does now 
take into account the oscillation of the data). There is 
also an additional bonus: the acquires statistical sig¬ 
nificance and indeed its distribution agrees very well 
with that expected for the function with one degree 
of freedom. Two variants are possible: (a) deriving Az 
from the difference between the estimated y(Z) and 
the measured T(Z), or (b) directly using the fourth 
parameter of the fit Az = do- We will call these two 
variants Y-COSa and Y-COSb, respectively. The re¬ 
sult of just one fit in the five-point interval around the 
charge Z=7 is shown by the thin solid line in fig.[Uc). 
The thick dashed line shows the parabolic part y{Z) 
in eq. ©• The estimate oi Az obtained with variant 
(a) is indicated once more by the arrow pointing from 
the experimental Y{Z) (open cross) to y{Z) (full cir¬ 
cle). It is apparent that now the circle represents a 
much better estimate of Ysm{Z) than in the previous 
method. 


5. Log-Third difference method (LY-3DI). This is actually 
the oldest method. First proposed by Tracy [2], it has 
been widely used in the past and also in recent pub¬ 
lications (see e.g. mm)- It relies on finite difference 
calculus (see Appendix A) and in particular it uses the 
third difference of the natural logarithm of the mea¬ 
sured yield L{Z) = \n{Y{Z)) over four consecutive Z 
values: 

4+1/2 = {L(Z + 2)-3L(Z+l) 

+3L(Z)- L(Z-l)}. (7) 

The basic assumption is that, without staggering ef¬ 
fects, the shape of the logarithm of the yield, L{Z), 
over the considered interval of four Z values could be 
described by a parabola 

L^Z) = ot Z'^ -\- bq- Z -\- ct- (8) 


The justification by Tracy was that the smooth be¬ 
havior of Ysni(Z) over a small interval of Z can be well 
described by a piece of a Gaussian. It is worth not¬ 
ing that this does not necessarily imply that the whole 
charge distribution needs to be Gaussian-shaped, as 
it was usually assumed in old-time low-energy fission 
studiesEI. One can easily verify that with the functional 
form of eq. ([5]) (and with no staggering effects) the fol¬ 
lowing identity holds 


ns 


D 


0) = Q 

Z+1/2 — " 


(9) 


^ All non-linear fits were performed by means of Minuit 
| 40 | (with successive calls to the routines Migrad, Hesse and 
Minos) and always reached convergence. 

^ Actually, as either sign is allowed for the coefficient ar in 
eq. (O, the shape of Y(Z) is not always a piece of a Gaussian. 



Fig. 2. (color online) Illustration of the methods (a) LY-3DI 
and (b) LY-4DI. Open squares and crosses are the same as in 
fig- [II The effect of the estimated staggering is shown by the 
arrows pointing from the experimental data (crosses) to the 
reconstructed smooth yields (open circles). The full lines are 
the assumed curves that best fit the open circles and the full 
circle gives the so estimated value of Ysm: see text for details. 
Note the expanded vertical scale with suppressed origin. 


If now a fixed quantity Dz is alternatively added to 
L{Z) for even Z and subtracted for odd ones, then 

one obtains T)^z+i /2 = easily recognize 

Wahl’s parametrization m of the staggering with Dz = 
hY{Fz), see eq. ([3]). The factor (—1)'^ in eq. 0 takes 
care of the fact that for staggering the sign of Dz is 
usually positive for even Z and negative for odd ones. 
We note also that the result is usually associated with 
the center of the interval and therefore, being D^^'> 
an odd finite difference, that center corresponds to a 
half-integer Z value. An illustration of this method is 
shown in fig.|2](a), where the arrows pointing from four 
consecutive experimental values of Y(Z) (crosses) to 
the corresponding open circles represent the estimated 
staggering effect over that interval and the solid line 
is the quadratic polynomial, in the logarithm of the 
yields L{Z) = \n(Y{Z)), that best fits the four open 
circles. The central full circle (red online) gives the so 
estimated value of at the center of the considered 
interval. 

The method LY-3DI by Tracy, applied to four consecutive 

points, is just a particular application of the mathemati- 
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cal formalism of Finite Difference Calculus, as briefly il¬ 
lustrated in Appendix A. This formalism has been widely 
used and discussed in the literature in connection with 
nuclear structure features, like pairing effects in nuclear 
masses and binding energies (see e.g. 131133), but not so 
much in connection with the fragment production yields. 

Actually, already the calculation for the first method, 
LIN, produces a 3-point equation that corresponds to the 
second-order difference in the fragment yields. For this 
reason we will henceforth address that method as Y-2DI: 
it represents the simplest useful implementation of this 
formalism (in fact, a first-difference 2-point formula would 
probably be too rough). 

Other implementations are possible, but - to our knowl¬ 
edge - they were not considered for studying odd-even ef¬ 
fects, at least not for what concerns the yields of reaction 
products. Applying eq. (jA4l) of Appendix A to the yields 
V(Z) and eq. (IA7I) to their logarithms L{Z), one obtains 
additional methods that are briefly sketched hereafter. 

6. Log-Second difference method (LY-2DI). For sake of 
completeness, we mention that a linear method, similar 
to Y-2DI, could be applied also to the logarithm of the 
yields, leading to the second-difference expression 

Df = (-1)^ \ {-L{Z-1) + 2L{Z) - L{Z+1)} 

( 10 ) 

which is again an estimator oi Dz, with the superscript 
indicating the order of the finite difference. As in the 
case of Y-2DI, also this method uses three points. 

7. Log-Fourth difference method (LY-4-DI). Much in line 
with the method of Tracy [5] , in the present paper we 
consider of great interest also the fourth difference of 
L{Z) over five consecutive points: 

i-lf ±{LiZ-2) -AL{Z-1) 

+ 6L{Z) -4L{Z + l) +L{Z+2)}. (11) 

Again the assumption is that, without staggering, L{Z) 
is well described by eq. ([8]) and also in this case one 
obtains = 0. However, if a quantity Dz is alter¬ 

natively added to and subtracted from L{Z), one finds 
again = Dz- With respect to the method if Tracy, 
this one uses five points (like all previous methods ex¬ 
cept Y-2DI and LY-3DI) and the staggering parameter 
Dz can be more physically attributed to the integer Z 
value at the center of the considered interval. An ad¬ 
ditional advantage of LY-4DI over LY-3DI is that eq. 
m remains valid even if L{Z) has some cubic com¬ 
ponent (see Appendix A): 

L{Z) = a Z^+ b Z"^+ c Z+ d. (12) 

An illustration of LY-4DI is shown in fig. HKb), with 
the same meaning of symbols and lines as in fig. Ha). 

8. Third difference method (Y-3DI). The difference of this 
method with respect to LY-3DI is that it uses the frag¬ 
ment yields (instead of their logarithms) of four frag¬ 
ments, giving a value of Az that is again attributed 


to the half-integer intermediate value Z From eq. 
dsn the staggering parameter A^ is obtained as 

^zll/2 = (- 1 )^ • 

{Y{Z-h2)-3Y(Z-hl)-h3Y{Z)- Y{Z-1)}/ 
{Y{Z+2) + 3 Y{Z + 1) + 3 Y{Z) + Y{Z-1)} (13) 

9. Fourth difference method (Y-4DI). The next imple¬ 
mentation with five fragment yields, again based on 
the fourth finite difference of eq. dSl), reads: 

= {-if {Y(Z-2) - 4 Y(Z-l) 

-h6 Y(Z) - 4 Y{ZF\) Y Y(Z + 2)}/ 
{Y{Z-2)y4Y{Z-\) 
+6Y{Z)+4Y{Z+l) + Y{Z + 2)} (14) 

This method has been extensively employed in nuclear 
structure studies of odd-even mass staggering m- 

Some comments are due, before beginning to examine 
the merits and drawbacks of the various methods that can 
be used to extract the staggering parameter {or Dz)- 
First of all, it is not advisable to use still higher-order finite 
differences because they perform weighted averages over 
an increasing range of experimental points. One has to find 
a compromise between using more points to be less sen¬ 
sitive to higher derivatives of the smooth yield function, 
and using less points to be more sensitive to local varia¬ 
tions of the staggering phenomena. We decided not to go 
beyond five points, which is the number commonly used 
also by the other methods based on fitting procedures. 

We also point out that the methods based on fits and 
those using finite differences of the fragment yields adopt 
the staggering parametrization by Amiel-Feldstein [I1II3], 
while those using finite differences of the yield logarithms 
find more natural to adopt that by Wahl [14] . We will see 
that this has only minor consequences. 

A further dissimilarity is that the methods relying on 
fits do first try to extract, from the data themselves, the 
smooth yield y{Z) (full circles) that best estimates the 
unknown Ysni(Z)- This in turn allows one to determine 
Az from the comparison with the experimental yield: 

Az = {-\)^ {Y{z)-y{z))/y{z) (15) 

(actually also Y-COSb gives directly from the fit). In 
eq. (ITKI) the factor (—1)^ takes care again of the opposite 
signs for even and odd Z values. 

On the contrary, in the methods using the finite dif¬ 
ference formalism there is generally no need to compute 
y{Z), because they give directly Az (or Dz), which is 
the quantity of physical interest. In figs. H^-) (b) the 

underlying Ys^{Z) distributions have been reconstructed 
just for illustrative purposes: by subtracting the stagger¬ 
ing effect from the experimental yields Y {Z) (crosses), one 
can deduce the smooth yields (open circles) and draw the 
curves (solid lines) corresponding to eqs. (|S|) and (jUj). 
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The two parametrizations by Amiel-Feldstein and Wahl 
are approximately equivalent in case of small staggering ef¬ 
fects (Az 1), being to first order Dz ~ In fact, tak¬ 
ing the ratio between the yields for even and odd species 
in both parametrizations, one obtains the relationship m 

(1 + Az)/{1 - Az) = (Fz) 2 (16) 

and expanding 1-1-2 -|- O {A\) = 1-1-2 Dz + O (-D|). 

In case of larger staggering, a fair quantitative compar¬ 
ison of the results obtained with the two parametrizations 
may require corrections of higher order m-- 

Dz = Az + I Az^ + i Az^ + O (Az'^) (17) 

Az = Dz-IDz^ + ^Dz^ + 0(Dz^) (18) 

3 Simulation 

A comparison of the different methods, to find out their 
possible advantages and disadvantages, is best performed 
by means of simulations. In this way one can play with 
the smooth yield distribution YsmiZ) and with the pa¬ 
rameter Az, order to produce a final distribution y(Z) 
(including statistical fluctuations) that simulates the ex¬ 
perimental one. This allows one to test in how far the 
different methods are able to retrieve the original Vsm(Z) 
and the genuine odd-even effect Az - In the simulations for 
the present paper, the Amiel and Feldstein parametriza- 
tion (Az) has been adopted and the results of the LY-2DI, 
LY-3DI and LY-4DI methods (that deliver Dz) have been 
transformed by means of eq. (HID. 

The simulation proceeds in the following steps: 

1. A smooth analytic function f(x) is assumed and its 
evaluation for integer x (x € N) gives the discretized 
smooth charge distribution Ysm{Z) without staggering. 

2. Following the parametrization by Amiel and Feldstein, 
the charge distribution with staggering is obtained as 
Y(Z) = YsmiZ) (1 -I- (-1)^Z\^), where in principle Ziz 
could even be a slow-varying function of Z. 

3. A proper normalization converts the yield Y(Z) into 
the number of counts Nc(Z) and then fluctuations 
are added to Nc(Z), thus producing an experimental 
charge distribution that simulates the outcome of an 
experiment. Statistical fluctuations are produced by 
means of the random number generator Ranlux m- 
In principle they should be Poisson-distributed, but 
the normalization is chosen to give values of Nc(Z) 
large enough, so that the fluctuations can be extracted 
from a Gaussian distribution with cF(Z) = Nc(Z). 

4. All the described smoothing methods are applied to 
the experimental distribution to obtain values of zl^ 
as a function of Z . 

Steps [3] and S] are repeated for a preset number of times, 
which in our case is always 10^. 

5. The propagation of the fluctuations on the individ¬ 
ual yields, including all possible correlations and non- 
linearities of the methods, produces - for each Z - 


a distribution of Az that can be well fitted with a 
Gaussian. The centroids (Az(Z)) give the distribution 
of the average result (over 10"^ distributions) for each 
method. 

6. The results obtained with the different methods are 
finally plotted and examined in detail. 

With reference to the above step [SI we emphasize that our 
procedure is equivalent (as we have easily verified with 
the simulation) to adding up all 10'^ replicas and applying 
the various methods to the summed distribution (which 
has a huge statistic). The advantage of our procedure is 
that the widths of the Gaussians, a a, give the uncertainty 
on Z\z - at the Icr level - expected in the analysis of 
one single experimental distribution. This is exactly the 
meaning that will be given to the dashed lines in the panels 
(b) to (h) of fig. island to the error bars in the panels (b) 
to (m) of figs, m to HI 

Our final aim is to compare the behavior of the dif¬ 
ferent methods of sect. HI in relation to both magnitude 
of the staggering effect and shape of the charge distribu¬ 
tion, as well as to assess the robustness of the result with 
respect to the statistics of the collected data. 

For what concerns the last point, to avoid that intrinsic 
differences among the various analysis methods are ob¬ 
scured by fluctuations, the relative errors on the yields 
should be kept appropriately small. In the simulation, 
where the fluctuations are just of statistical nature, it 
is enoi^h to use a reasonably large number of counts 

NciZ)E 

One should be aware that the fluctuations on one Z 
affect the estimated value of the staggering parameter Zl^ 
also for neighboring charges, and this happens in quite a 
complicated and correlated way that is difficult to esti¬ 
mate without a simulation. In fact, while in a real exper¬ 
iment one can analyze just one single realization of the 
charge distribution, a big advantage of the simulation is 
that one can generate a large number of replicas of the 
same charge distribution (differing from each other only 
in the random fluctuations) and all correlations come nat¬ 
urally out from the simulation itself. 

Finally, the simulation allows one to choose freely the 
shape of YsmiZ) in order to see the response of the differ¬ 
ent methods to the shape of the charge distribution, with 
special attention to its regions of non-linearity. 


4 Results 

The first thing that we checked is that indeed all meth¬ 
ods are unbiased, i.e. in absence of staggering they give 
Az = 0 within errors. This is shown in fig. [31 where we 
consider constant charge distributions, without stagger¬ 
ing, with a statistic of 200, 500 and 1000 counts per Z 
unit (first, second and third columns, respectively). In the 

^ In a real experiment this may not suffice, as one has to in¬ 
clude other possible uncertainties, if any, that may affect neigh¬ 
boring charges Z in an independent way (e.g. contaminations, 
identification uncertainties, and so on). 
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Fig. 3. (color online) Effect of statistical fluctuations for three 
constant charge distributions with 200, 500 and 1000 counts 
per Z unit (first, second and third column, respectively), with¬ 
out staggering {Az = 0). Panels in first row (a): relative sta¬ 
tistical fluctuations, (F- (Y'))/(y), generated by three differ¬ 
ent sequences of random numbers (full, red dashed and dotted 
lines), that applied to the constant distributions produce three 
experimental realizations. Panels in rows from (b) to (h): av¬ 
erage values (full dots) of the parameter Az evaluated over 
10"^ distributions with some of the methods illustrated in this 
paper, compared with Az obtained for each of the three real¬ 
izations (see text). 

first row, fig. ED a), the continuous, dashed (red online) and 
dotted lines show the relative errors AYjY obtained by 
applying the same three sequences of random errors (step 
3 of sect. ED to the above mentioned constant distribu¬ 
tions: in fact the structure is the same, the relative errors 
being reduced just because of the increasing number of 
counts. The horizontal dashed lines show the ilcr error 
band expected for just one realization and correspond to 
relative errors AYjY of about ± 7% for Nc=200, down to 
about ± 3% for iVc=1000. This means that in one single 
realization of the charge distribution (i.e. in one measure¬ 
ment), for each measured Z there is still a probability of 
~ 32% that its measured yield lies outside this band. 

These relative errors may seem rather small, but in 
the evaluation of Az they are amplified by all examined 
methods. In fact propagating these errors from the yields 
to Az leads to the ±lcr band indicated by the horizontal 


dashed lines in rows (b) to (h) of fig. El its width amount¬ 
ing for all methods to about ± 0.04 when Nc=200 and ± 
0.02 when Afc=1000B 

But this is not the worst part of the story. The impact 
of these fluctuations on the estimation of Az is presented 
in the lower panels of fig. El rows from (b) to (h). One 
clearly sees that taking the average (step 5 of sect. ED over 
10"'^ generated replicas, differing only in the fluctuations, 
gives = 0 with all methods (full dots). However, the 
correlated propagation of the fluctuations of neighboring 
charges gives rise to rather large and nonphysical struc¬ 
tures (full, dashed and dotted lines), with different and 
often opposite phases in different realizations: one realiza¬ 
tion may have a bump where another one has a valley. It 
is also worth noting that the spurious bumps (or valleys) 
depend only on the particular set of random fluctuation 
and appear to be strictly in phase for all analysis methods 
that can be applied. 

The main conclusions that can be drawn from fig. El 
are that (a) fluctuations in the charge distribution are 
amplified by all methods and the values of Az for neigh¬ 
boring Z become strongly correlated; (b) good statistic 
(and reduction of other sources of fluctuations) is there¬ 
fore mandatory for any experiment addressing this topic; 
(c) structures or trends in Az{Z) should be interpreted 
with great caution, after a careful and realistic estimate 
of their statistical significance. 

In the remaining of this paper we will adopt a constant 
and representative value Az = 0-1. In fact we played with 
possible variations of the staggering magnitude and of its 
charge dependence, but nothing noteworthy was found. 
On average (or, equivalently, in a single realization with 
very small relative errors) almost all methods do an ex¬ 
cellent job when the staggering Z\z = 0.1 is applied to a 
charge distribution that either stays constant, or rises or 
falls in a linear way, whatever the value of its slope. The 
only exception concerns the method Y-PAR that, with 
respect to the other methods, amplifies the staggering by 
about 40% (see Appendix B for more details). 

Apart the trivial (and unrealistic) case of a linear charge 
distribution, we now have to investigate how the different 
methods behave in regions where the distribution Ys^(Z) 
is highly not linear (namely it presents an appreciable cur¬ 
vature or has an inflection point) on the basis of some 
representative cases that are shown in figs. E] to [71 

The figures are organized as follows. The first panel (a) 
always shows the yield distribution, while the remaining 
panels, from (b) to (m), present the result (Az) obtained 
with the various methods and averaged over 10'^ repli¬ 
cas of the charge distribution (from now on we will omit 
the ( ) for simplicity). The results in the first three rows 
(panels from (b) to (g)) are based on the second, third 
and fourth finite differences, applied either directly to the 
yields (Y-2DI, Y-3DI and Y-4DI in the left columns), or 
to their logarithms (LY-2DI, LY-3DI and LY-4DI in the 
right columns). The results in the last two rows (panels 

® Similar absolute uncertainties occur also in presence of 
staggering and cannot be overlooked: for Az =0.1 they corre¬ 
spond to an uncertainty of « ±40% (20%) for Ac=200 (1000). 
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Fig. 4. (color online) Upper panel (a): charge distributions shaped as vertical parabolas (pairwise with the same curvature). 
Open dots show a constant staggering Z\z=0.1 applied to the smooth dashed curves. Panels (b)-(m) present the average 
{Az) obtained with the methods Y-2DI, LY-2DI, Y-3DI, LY-3DI, Y-4DI, LY-4DI, Y-2DIF, Y-PAR, Y-COSa and Y-COSb, 
respectively, compared with the nominal input value (dotted lines). The error bars are the Tier widths of the Az distributions. 


from (h) to (m)) always refer to the methods based on 
various types of fits to the yield distribution (Y-2DIF, Y- 
PAR, Y-COSa and Y-COSb). Note that the vertical scales 
for the methods Y-2DI and LY-2DI (panels (b) and (c)) 
and for Y-2DIF and Y-PAR (panels (h) and (i)) span a 
range about twice as large as that for the other meth¬ 
ods. As already anticipated, the error bars represent the 
statistical uncertainty (at the ilcr level) for the analysis 
of just one realization of the distributions shown in the 
panels (a). 

Figure m shows six examples (three maxima and three 
minima) of charge distributions shaped as vertical parabo¬ 
las y{x) = a(x — xq)^ + c. They have pairwise the same 
curvature, but with opposite sign (a = ±200 for the first 
pair, ± 400 for the second and ± 600 for the third one), 


while the value in the vertex is always the same, c = 20000. 
With a staggering parameter of Az=0.1, the enhancement 
(reduction) of even (odd) charges Z near the vertex is al¬ 
ways the same, i.e. ±(-)2000 counts (step 2 in sect.[31). 

As expected, the simple Y-2DI and LY-2DI, shown 
in panels (b) and (c), are not well suited for distribu¬ 
tions with an appreciable curvature and indeed they give 
the worst results. For what concerns the remaining finite- 
difference methods they all behave rather well. The method 
LY-4DI in (g) gives some improvement over LY-3DI, the 
method by Tracy shown in (e), especially for the down¬ 
ward parabolas at the two extremes of the usable Z range 
(where the statistic is very low). However the linear meth¬ 
ods Y-3DI (d) and Y-4DI (f) seem to be slightly better. 
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Fig. 5. (color online) Similar to fig. [d] but for charge distributions with different shapes of the inflection points. 


For the methods based on fits to the data, the ad¬ 
dition of two more points to Y-2DI that is necessary to 
perform a fit (method Y-2DIF in (h)) certainly improves 
the results, especially for the smallest curvatures. In the 
Y-PAR method, panel (i), one sees the already mentioned 
amplification of about 40% and the intrinsic difference of 
even and odd Z (already reduced by treating the errors 
as explained in Appendix B). Y-PAR gives similar results 
for all curvatures (as expected from its being based on a 
parabolic fit), but it certainly does not represent a real im¬ 
provement over Y-2DIF. Much better results are definitely 
obtained by the two new methods, Y-COSa and Y-COSb 
(based on the ht with a properly oscillating function). 

The situation is in general even worse for inflection 
points, which make more trouble for nearly all methods. 
In fig. [SKa) we present six charge distributions that differ 
either in the slope at the inflection point or in the distance 


between the minimum and the maximum. The distribu¬ 
tions are obtained from cubic functions without quadratic 
term, Y(Z) = a (Z—Zo)^+b (Z—Zo)+ c, all with the same 
statistic at the inflection point (c = 20000). The slope at 
the inflection point is 6 = 1000, 2000 and 3000 for the first, 
second and third pair of distributions, respectively, while 
the cubic coefficient is a = —50 for the first and —150 for 
the second curve of each pair. 

The lower panels of fig-El from (b) to (m), present the 
average results for the different methods. The methods Y- 
2DI and LY-2DI of panels (b) and (c) are certainly not 
acceptable, because they give systematic deviations up to 
±0.05 with respect to the nominal value of Az=0.1. Of the 
remaining finite-difference methods, Y-4DI and LY-4DI of 
panels (f) and (g) are now very remarkable improvements 
over Y-3DI (d) and LY-3DI by Tracy (e), because they 
reduce the deviations from the nominal value down to less 
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Fig. 6. (color online) Upper panel (a): high-statistic charge distribntion based on a natural cubic spline (smooth dashed line); 
the thin vertical lines indicate the internal knots (see text). Open dots display the effect of a constant staggering Az = 0.1. 
Panels (b)-(m) present Az, as obtained with the methods discussed in this paper, compared with the nominal value of 0.1 
(dotted lines). For each Z, the error bar indicates the ilcr width of the Az distribution obtained from 10"^ simulations. 


than ±0.005. This does not come totally unexpected, be¬ 
cause methods based on the fourth differences should be 
able to remove smooth terms up to the third order (see 
Appendix A). 

The fit method Y-2DIF of panel (h) represents again 
an appreciable improvement over the linear method Y-2DI 
(b), as it cuts the systematic deviations down by about a 
factor of two, but the method Y-PAR of panel (i), in spite 
of its amplification and systematic even-odd effect, per¬ 
forms even better than Y-2DIF. What comes rather un¬ 
expected is that the methods Y-COSa (1) and especially 
Y-COSb (m) perform rather well, in a way comparable 
with Y-4DI (f) and LY-4DI (g). Altogether, the best re¬ 
sults among all those displayed in fig. [S] are obtained with 
the methods Y-4DI and Y-COSb. 


The high statistic adopted for the smooth distribution 
Ysm{Z) (with relative errors lower than ±1% in the non¬ 
linear regions) simply makes the presentation of the data 
in fig. m and [S] more readable by reducing the error bars, 
but obviously does not appreciably change the response 
of the methods, averaged over many distributions [j. Of 
course, for a single realization (or experiment) a large 
statistic of the charge distribution remains of paramount 
importance. 


® For example, one sees from eqs. (HU, dl) or m that 
a common factor k multiplying all yields cancels out in the 
average {Az)', and it cancels out also in (Dz), when (ln(fc Y)) = 
In(fc) + (ln(y)) is inserted in eqs. (0 or m- 
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Fig. 7. (color online) Upper panel (a): same charge distribution as in fig.[6l but with low statistic. In panels (b)-(m) the average 
behavior over 10^ realizations (full dots) is the same as in the previous figure, only the fluctuations (error bars) expected for 
a single realization are much larger. The open dots/squares (joined by dashed/dotted lines) put into evidence the correlated 
propagation of the statistical fluctuations for two randomly selected realizations of the low-statistic charge distribution. 


A second point worth noting is that the deviations 
from the nominal value displayed in fig. |4] and [5] are char¬ 
acteristic of each method and of the specific non-linearities 
of the smooth function The simulation also shows 
that in general the deviations from the nominal value of 
the staggering parameter are very weakly dependent on 
the chosen value of For methods based on the yield 
logarithms, one can even demonstrate that these devia¬ 
tions do not depend on it. Considering for example the 
case of LY-3DI, if L{Z) has the form of eq. (TT^ with a 
cubic coefficient a, then eq. © gives Dz + fa (instead 
of Dz), thus with a constant deviation from the nominal 
value that amounts to |a and hence does not depend on 
Dz- Therefore, the choice of a good method is important, 
especially when the average staggering effects are small. 


We tested several other shapes of the charge distribu¬ 
tion, finding that in all cases, irrespective of the specific 
shape, the ranking of the various methods remains similar 
to that found in figs. [Hand[S] 

We show just one case in figs. H] and [3 where a more 
realistic charge distribution - open dots in panels (a) - 
mimics the typical features found in semi-peripheral colli¬ 
sions of intermediate-energy heavy-ion reactions: with in¬ 
creasing Z there is a very rapid drop of the cross section 
from the lightest complex fragments with Z ks 4 down to a 
wide valley of intermediate mass fragments (IMF) around 
Z Ri 10 - 15, followed by a slow rise towards heavier ones 
(which may represent either fission-like or projectile-like 
fragments) until it finally fades away. 
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This shape was obtained with a natural cubic spline 
passing through seven equally spaced knots (the five inter¬ 
nal ones are indicated by thin vertical lines in the figures): 
it is a piece-wise interpolation with 3'’*^ degree polynomials 
that is continuous in the joining knots up to the second 
derivative. The adjective natural indicates the additional 
condition that the second derivative is zero in the two 
external knots. 

In fig.[6]the statistic for each of the 10"* generated repli¬ 
cas is very large: in the IMF valley it corresponds to at 
least ten thousand counts for each charge Z. The lower 
panels, from (b) to (m), illustrate again the estimates of 
Az obtained with the various methods already shown in 
the previous figures. The full dots represent the average 
obtained over 10^ simulations, while the error bars, often 
smaller than the symbol sizes, give the widths of the dis¬ 
tributions for each Z. Again dashed lines at Az = 0.1 
indicate the nominal value. 

Also here we observe an evident failure of Y-2DI, LY- 
2DI and Y-2DIF in the minimum of the IMF valley and 
in the final tail of the distribution. Of all finite-difference 
methods, those based on the fourth differences perform 
better (Y-4DI even slightly better than LY-4DI), but com¬ 
parably good results are obtained with the fitting methods 
Y-COS and particularly with the second one, Y-COSb. 

Concerning the finite difference methods, one might 
think that fragment cross sections that vary over many 
orders of magnitude are better studied with the logarith¬ 
mic n-point differences. The present study does not lend 
much support to this idea, as the results of linear and log¬ 
arithmic methods presented in figures from 0] to [B] do not 
differ too much from each other. For the second-difference 
results of panels (b) and (c), one might have a slight pref¬ 
erence for the logarithmic methods, but the choice is hard. 
The situation for the third and fourth differences (panels 
(h) vs. (i) and (i) vs. (m), respectively) seems to be slightly 
in favour of the methods that use directly the yields. In 
fact, the logarithmic methods seem to suffer more in re¬ 
gions where the yield is lower, see e.g. the IMF valley of 
fig. m or the sides of the downward parabolas in fig. 01 

In order to stress the importance and the effects of the 
statistic, fig. [7] presents results for the same charge distri¬ 
bution of fig. 01 but in the case of low statistic: about 300 
counts per charge Z in the IMF valley. The average be¬ 
havior is substantially the same as in the previous figure. 
In fact, for each of the presented methods, the full dots are 
practically the same in both figures, their differences being 
much smaller than the symbol sizes. The main difference 
concerns the ±lcr fluctuations (error bars) expected for a 
single realization (or measurement), which are of course 
much larger in fig. [T] 

In addition, two examples of the results obtained for 
two individual random realizations of the charge distri¬ 
bution are shown by the open dots/squares, joined by 
dashed/dotted lines in the panels (b)-(m) of fig. [T] Here 
the effects of the statistical fluctuations are clearly vis¬ 
ible for all methods and reveal marked deviations from 
the average behavior, with several points located outside 
the error bars. What is even worse is the fact that the 


fluctuations give rise to structures involving several neigh¬ 
boring points, a consequence of their correlation (see next 
section). The location and the amplitude of these non¬ 
physical structures present large random variations from 
one single realization to the other, being sometimes in- 
phase and sometimes out-of-phase in the two realizations. 
This means that in a single experiment, unless the inde¬ 
pendent relative errors, including the statistical ones, are 
very small (below the 1% level), the physical significance 
of structures in the staggering parameter Az should be 
treated with great caution. 


5 Correlations 

The need of estimating the smooth behavior of the charge 
distribution from a few data points, naturally introduces 
short-range autocorrelations among the staggering param¬ 
eters Az pertaining to neighboring charges Z. The degree 
of autocorrelation is estimated with the usual definition 

R{Z', Z") = - Az-')) 

(TAz, OAz„ 

where the overlined Zl^’s are the mean values of the stag¬ 
gering parameters for two different charges Z' and Z" ^ the 
brackets ( ) indicate the expectation value of the product 
of the deviations of the two Az from their mean values and 
the cr’s in the denominator are their standard deviations. 
By construction R{Z', Z") lies in the range [—1,1], the two 
extremes indicating perfect correlation or anticorrelation, 
respectively (while 0 indicates uncorrelated data). 

The correlation factor R cannot be estimated from a 
single experiment, but it is easily obtained in a simulation 
where many replicas of a certain charge distribution are 
generated, differing only in the fluctuations. As expected, 
fluctuations introduce a correlation that is positive for two 
nearby located Z' and Z" values and 0 for distant ones. 

The correlation factor is practically independent of the 
specific shape of the analyzed charge distributions. For a 
flat distribution with a constant Az of 0.1, the region 
of positively correlated neighboring Z values is obviously 
at most three charge units wide for the methods Y-2DI 
and LY-2DI (with FWHM values of 2.37±0.03) and in¬ 
creases to four charge units for Y-3DI and LY-3DI (with 
FWHM values of about 2.84) and to five charge units for 
all remaining methods (with a common FWHM of 3.24). 
In conclusion, structures possibly found in the staggering 
analysis that are narrower than the above numbers are 
highly suspicious and probably they should not be given 
much physical significance. 


6 Conclusions 

The odd-even staggering is a widespread phenomenon that 
usually consists in a relative enhancement of the yields of 
nuclei with even values of proton (neutron) number Z (N) 
with respect to the neighboring nuclei with odd values. All 
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methods are based on the idea, first proposed by Amiel 
and Feldstein that the measured yield Y{Z) can 

be factorized into the product of a smooth yield Ys^{Z) 
times a factor (1 ± Z\z), where the sign ‘-I-’ holds for even 
charges and the sign ’ for odd ones. 

This is an assumption, widely accepted and probably 
justifiable on the basis of theoretical arguments, but it 
remains an assumption and one should definitely remain 
aware of that. 

We have reviewed some old methods and we propose 
some new methods of extracting the staggering parameter 
Az (or Dz) from the experimental charge distributions. 
Some methods first estimate the underlying smooth be¬ 
havior of the yield y{Z) from the data themselves, to de¬ 
rive then the staggering from the difference between the 
measurement and the smoothed yield, while others de¬ 
liver a direct estimate of the staggering parameter. Some 
methods use a fit to the data, others take advantage of 
the finite-difference mathematical formalism. 

Almost all methods perform equally well if the charge 
distribution has a predominant linear behavior, while they 
greatly differ from each other in the presence of strongly 
non-linear regions in the charge distribution. So steep¬ 
ness is not an issue, not even for fragment yields. From 
the present study, using logarithmic n-point differences 
for fragment yields is just an option, not a real necessity. 

All methods are extremely sensitive to the fluctuations 
that may affect the yields of neighboring Z in the experi¬ 
mental data, giving origin to spurious structures without 
much physical meaning. Therefore for this kind of studies 
it is of paramount importance to acquire good data with 
very large statistic and very small independent relative 
errors, so that the relative random fluctuations for each 
measured Z are of the order of about 1% or better: only 
when this condition is fulfilled, one can really appreciate 
the different quality of the various methods. On the con¬ 
trary, if the fluctuations are sizably larger, spurious effects 
appear, and it makes no much sense to choose one method 
rather than another; the results will be anyhow poor and 
hard to interpret in a sensible way. 

We find that linear interpolations (Y-2DI and LY-2DI) 
that make use of three neighboring points give the worst 
results among all examined methods. Extending the lin¬ 
ear interpolation with a fit that uses five points (Y-2DIF) 
brings only some moderate improvement. The recently 
proposed fit of a simple parabola over five points (Y-PAR) 
does not represent a real improvement over Y-2DIF. The 
method Y-PAR presents the peculiarity that it amplifies 
the signal, but it needs a careful treatment of the errors 
and even then it does not produce a really smooth esti¬ 
mate of the distribution without staggering, y{Z). 

Very good results are given by the newly proposed 
method (Y-COS) that uses a five-point fit with a prop¬ 
erly oscillating function. As it implies a non-linear fit, a 
reliable recursive procedure has to be used, such as the 
well known Minuit code [JD] developed in CERN. Using 
directly the coefficient of the oscillating term (method Y- 
COSb) as delivered by the fit seems to be one of the best 
ways to estimate the staggering parameter Az- 


A completely different approach is the oldest one by 
Tracy, based on the third (LY-3DI) finite differences over 
four values that are the logarithms of the fragment yields. 
At variance with the previous methods, it delivers directly 
an estimate of the staggering parameter, although with a 
slightly different parametrization first proposed by Wahl 
|14) . We show that a similar method (Y-3DI) could be ap¬ 
plied to the fragment yields, using the parametrization by 
Amiel [HIS]. In both cases the staggering value, that is 
usually attributed to the center of the considered interval, 
corresponds to a half-integer Z. These two methods give 
good results, but not as good as those of Y-COSb. 

Eurther exploiting the Einite Difference formalism, in 
this paper we propose two more methods, Y-4DI and LY- 
4DI, that use the fourth finite differences over five points. 
They give somewhat better results than the original LY- 
3DI because they are insensitive to possible cubic com¬ 
ponents in the smooth charge distribution, and probably 
also because they use one more point and hence they are 
slightly less sensitive to fluctuations. Moreover, one prac¬ 
tical (and aesthetic) advantage is that the results can be 
attributed to the integer Z values at the center of the 
five-point intervals. 

Concluding, to avoid the appearance in Az{Z) of struc¬ 
tures without much physical significance i) one should 
collect data with a very large statistic; ii) all other un¬ 
certainties that may affect the yields of neighboring Z in 
an independent way should also be carefully estimated 
and reduced; iii) once the previous conditions are met, 
the methods which seem less sensitive to non-linearities of 
the underlying smooth distribution YgmiZ) are the meth¬ 
ods Y-4DI and Y-COSb, closely followed by LY-4DI and 
Y-COSa. 


Appendix A 


The Einite Difference Calculus (see e.g. mMM) is a 
widely used mathematical formalism dealing with finite 
increments of the independent variable(s) of mathemati¬ 
cal functions. In particular, here we consider a function 
y{x) whose values are known for a series of equidistant 
values of the independent variable x; such is the case of 
the fragment yield Y{Z) that is defined only for integer 
values of Z. 

The known values of the function form a sequence yk , 
where the subscript fc = 0,1, 2 ... indicates the number of 
equal increments with respect to a given reference value xq 
(in the general case k can take negative integer values too). 
Erom the sequence yk one can build the first difference 
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Vyk = {yk+i — Uk), and all differences of higher order Q 
V^yk = - V'^-^yk) = E(-l)* 2/fe+z 

i=0 

(Al) 

where (”) are the binomial coefficients. 

The property relevant for the present paper is the fol¬ 
lowing: if the sequence yk is described by a polynomial of 
degree n-1 in the index k, namely yk = Cik^), then 

the differences of order n (built from n-l-l values of yk) are 
identically zero (and so is for those of higher order too) 0- 
If one further assumes that the sequence yk, described 
by a polynomial of degree n-1, is perturbed by a staggering 
parameter A a la Amiel 

71 — 1 

J/fe = (^c,f)(i + (-1)M), (A2) 

then a little algebra gives 

n z V 

(A3) 

i=0 

and hence the (n-l-l)-point formula: 


^k+n/2 = (-1)" 'D’^Vk / (E 

(A4) 

where we have added to A the subscript k+n/2 indicating 
the central value of the used interval [yk, Vk+n] ■ 

Assuming that the polynomial of degree n describes 
the logarithm of the sequence yk, it results more convenient 
to use the staggering model a la Wahl and according to 
eq. ([2]) one can write 

n —1 

ln(j/fc) = 

(A5) 

which, with a little algebra, gives 


V^yk = (2^) ■ D 

(A6) 

and again an (n-l-l)-point formula for D: 


Dk+u/2 = (-1)'= ^ V-yk. 

(A7) 


where again the subscript k-\-n/2 indicates the central 
value of the used interval. 

^ In mathematical textbooks the difference and derivative 
operators of order n are usually represented with A” and I?", 
respectively. Unfortunately, the same symbols are commonly 
used for the staggering parameters of Amiel/Wahl. For the sake 
of clarity, here we use X>"' for the finite difference operator. 

® If yk is generated by a function that is not a polynomial, 
then the finite difference operator D" can be used to remove 
all smooth terms in the Taylor series up to the order n-1. 


Identifying the fragment yield Y{Z) with yk (and with 
a trivial index redefinition) these equations generate all 
the methods of Sect. [5] that do not make use of a fit to the 
yield data. For example, from eq. (EZl with n = 3, using 
the dehnition of eq. m, one obtains the method LY-3DI 
of eq. dS]), first proposed long time ago by Tracy [5]; and 
with n = 2 (4) one obtains the method LY-2DI (LY-4DI). 
Similarly, from eq. (IA41) with n = 2, 3 and 4, one obtains 
the methods Y-2DI, Y-3DI and Y-4DI, respectively. 


Appendix B 

With respect to the others, the method Y-PAR ampli¬ 
fies the estimated Az by about 40% above the nominal 
value. This has to be taken into account when quantita¬ 
tively comparing the results of this method with those of 
the others. A second characteristic intrinsic in the Y-PAR 
method is that this amplification is not the same for even 
and odd charges. It is somewhat larger for even charges 
(usually enhanced by staggering) than for odd ones, so 
that even in case of constant input values of the yield 
Ysni{Z) and of the parameter Az, one does not obtain 
a really flat Az distribution, even for an input value of 
Az=0 (see full dots in panels of row (d) in fig. [3]). 

This effect is connected with the weights used in the 
weighted least squares fit procedure. When the smoothed 



Fig. 8. (color online) Effect of applying a weighted least square 
fit in the Y-PAR (a) and Y-COS (b) methods, with the weights 
taken from the statistics of the data. The original distribution 
(squares) is flat with a statistic of 1000 counts; a Az=0.1 stag¬ 
gering effect (crosses) is superimposed on it. In the Y-PAR 
method the differing statistical errors for even and odd Z in¬ 
duce an artifact that estimates a larger staggering for even Z 
than for odd ones (see text). 



















Alessandro Olmi, Silvia Piantelli: Methods for a quantitative evaluation of odd-even staggering effects 


15 



I .... I .... I .... I .... I .... I .... I .... I .. I I 

5 10 15 20 25 30 35 40 ^ 


Fig. 9. (color online) Staggering parameter Az for the charge 
distribution of hg. [TJa) obtained with the method Y-PAR, ns- 
ing the weighted least square fit with statistical errors estimated 
directly from the experimental data (open squares), or by first 
applying the LIN/Y-2DI linear method (full dots). 


yield y{Z) is deduced by means of a fit (like in the Y- 
PAR and Y-COS methods), some care must be put in the 
treatment of the errors on the measured yields Y(Z). In 
fact the weighted least square fit method may introduce a 
subtle asymmetry between even and odd Z values when 
the weights mainly depend on the number of counts mea¬ 
sured for each Z . In this case, even Z yields (which are en¬ 
hanced) have systematically somewhat larger errors (and 
hence smaller weights) than odd Z yields. 

This causes an artifact in the Y-PAR method because 
the parabola is not the appropriate function to fit the five 
‘zigzagging’ experimental points. This is shown in fig. IHa), 
where a constant yield of 1000 counts with Az=0.1 is as¬ 
sumed, leading to statistical weights of about I/vTTOO 
and l/VOOO for the even and odd Z, respectively. With 
their larger weights, the three reduced yields (lower crosses) 
of an odd Z are more effective in shifting the parabola 
downward than the three enhanced yields (upper crosses) 
of an even Z in shifting it upward, due to their smaller 
weights: as a consequence, even Z values appear to have a 
larger staggering than odd ones. We stress that this effect 
cannot be cured by increasing the statistic of the measure¬ 
ment. It is worth noting that, on the contrary, the fits of 
the Y-COS method shown in fig. |SKb) are not affected by 
this problem, because this method uses a properly oscil¬ 
lating fit function that well reproduces the data. 

The open squares of fig. [9] show the effect on Az when 
the errors in Y-PAR are estimated directly from the data, 
in case of the charge distribution already presented in fig. 
[TKa). This artifact may be cured either by recurring to the 
unweighted least square method (which is however unsatis¬ 
factory when the yields strongly vary with Z) or perform¬ 
ing appropriate averages of the weights before fitting the 
data. In this paper we always used as statistical weights 
for the Y-PAR method the results of a preceding analy¬ 
sis performed with the Y-2DI linear method. This artifice 
proved essential to obtain results (full dots in fig. [21) that 
could be meaningfully compared with the other methods. 

We wish to warmly thank Prof. P.R. Maurenzig for 
continuous support, useful suggestions and very fruitful 
discussions. 
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